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ABSTRACT 


The intensity of turbulence and the Lagrangian correlation coefficient 
for a LOX-GH 2 rocket combustion chamber have been determined from experimental 
measurements of tracer gas diffusion. A combination of Taylor’s turbulent 
diffiision theory and a numerical method for solving the conservation equations 
of fluid mechanics was used to calculate these ‘cuantities. Taylor's theory 
was extended to consider the inhomogeneity of the turbulence field in the axial 
direction of the combustion chamber, and an exponential fi;, .ntion of the form 
v;as used to represent the Lagrangian correlation coefficient. 

The results Indicate that the value of the intensity of turbulence 
reaches a maximum of 14% at a location about 7” downstream from the injector. 
The Lagrangian correlation coefficient associated v;ith this value is given by 
the above exponential expression where a = 10,000 sec~^. 
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SUMMARY 


The intensity of turbulence and the Lagrangian correlation coefficient 
in a LOX-GH^ rocket combustion chamber have been determined experimentally 
using a small rocket engine operating at a nominal chamber pressure of 
150 psia. The experimental method consisted of injecting a tracer gas 
at a point along the chamber centerline while taking samples along a 
diameter at a downstream station. Three sample stations were investigated. 

For each of the sample stations, several tracer gas injection points were used. 

The gas samples were analyzed using a combination of weight analysis, Orsat 
analysis aiid gas chromatography to determine the tracer gas concentration 
profiles. The turbulence parameters were then calculated from the tracer 
gas concentration data using a combination of 6. I. Taylor's turbulent 
diffusion theory and a numerical procedure developed by Rocketdyne for 
solving the conservation equations of fluid mechanics. Taylor’s theory 
was extended to consider the inhomogeneity of the turbulence field in 
the axial direction of the combustion chamber. The Lagrangian correlation 
coefficient was represented by an exponential function of the form e“ > 
where x is the difference in dispersion times for an individual fluid particle 
and a is a constant. 

The results indicate that the value of the intensity of turbulence 
reaches a maximum of 14% at a location about 7 inches do\mstream from the 
injector. The Lagrangian correlation coefficient associated with this value is 
given by the above exponential form where a =; 10,000 sec“^. This represents a 
relatively small degree of correlation in that the correlation will drop to 
above 10% in one inch of chamber length. 
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I. INTRODUCTION 


Turbulent mixing is one of the most important physical processes 
occurring in a rocket combustion chamber. It influences the combustion 
process and therefore the performance of the rocket engine. In 
turbulent flow theory, the turbulent mixing of fuel and oxidizer is 
quantitatively expressed in terms of the intensity of turbulence and 
the Lagrangian correlation coefficient. At present there are very 
little data available regarding these quantities which can be used in the 
design and development of rocket engines. 

For the past six years, a series of experiments has been conducted at 

Tulane University to determine the intensities of turbulence and the 

Lagrangian correlation coefficients in rocket combustion chambers using 

different combinations of fuel and oxidizer. The method used consisted of 

experimentally measuring the diffusion of a tracer gas. The Intensity of 

turbulence and the Lagrangian correlation coefficient were calculated from 

the diffusion measurements using a combination of Taylor's turbulent 

diffusion theory and a numerical procedure for solving the conservation 

equations of fluid mechanics. In the calculations, the turbulence field was 

assumed to be isotropic throughout the combustion chamber and homogeneous 

only in the radial and circumferential directions but not in the axial 

direction. The Lagrangian correlation coefficient was assumed to be an 

exponential function of the dispersion time. In earlier work, the intensities 

of turbulence and the Lagrangian correlation coefficients in the combustion 

(2) fa) 

chambers of a small LOX-heptane rocket and a GOX-GH^ rocket '• vrere 
determined. The present v/ork is the latest in the series of the experiments 
conducted. 

The objective of the present work is to investigate the intensity of 
turbulence and the Lagrangian correlation coefficient in a liquid oxygen (LOX) 
and gaseous hydrogen (GH^) rocket combustion chamber. These quantities were 
. experimentally determined by a tracer gas diffusion method similar to that 
used in the earlier work. A combination of Taylor's turbulent diffusion 

( 4 ) 

theory and a numerical method for solving fhe conservation equations of 

fluid mechanics was used to calculate the intensity of turbulence and the 
Lagrangian correlation coefficient. 


II. DESCRIPTION OF EXPEPJtlENT 


(1) Experimental Apparatus 

A small rocket engine, as shown by Figure 1, using liquid o>Qrgen and 
gaseous hydrogen and operating at a normal chamber pressure of 150 psia was 
used for the experijnents . The propellants were injected in the form of co- 
axial jets from seven ports in the injec*?.;:«r. Helium was used as the tracer 
gas and could be injected at various points along the centerline of the 
combustion chamber through the helium injection strut. A sample probe passing 
across the combustion chamber along a diameter v;as used to withdrav/ gas samples. 
This probe contained six sample ports and could be moved laterally across the 
chamber or rotated circumferentially about the chamber centerline. Further, 
the sample probe could be placed at several longitudinal stations. 

The combustion chamber walls, the tracer gas injector and the sample 
probe are water-cooled. 

(2 ) Experimental Procedure 

An automatic timing device was used to fire the rocket engine and take 
the samples. About 11/2 seconds are required for the engine to reach 
equilibrium. During this time the sample lines were vented to the atmosphere 
and the helium tracer gas was injected so that equilibrium condition could be 
reached. Samples v/ere withdravrn for about one second and collected in sample 
bottles for later analysis. 

(3) Analysis of Samples 

A method combining weight analysis, absorption methods and gas chromato- 
graphy was used to analyze each sample containing water, oxygen, hydrogen and 
helium. More than 90% of the sample was water whose weight could not be deter- 
mined either by absorption or by gas chromatography. Therefore, the sample 
bottle vras weighed before the analysis, and again after the analysis had been 
completed and the bottle evacuated. This allowed the total weight of the 
sample to be determined. 

Before the analysis was to begin, the sample bottle was charged with argon 
to a pressure of 20 psig at room temperature. V/ith the pressure and temperature 
of the sample plus argon knov/n, the mass of each component could be calculated 
using the volume fraction information from the gas analysis described belovf. 
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The volume fraction of ox.ygen ai\d most of hydrogen v/ere measured by 
selective absorption using an Orsat ana.yzer. The absorption method was used 
because it was well suited for analyzing these gases in large concentrations. 

The remaining sample which contained argon, helium and a small amount of 
hydrogen was analyzed using gas chromatography to determine the concentx-ation 
of helium and hydrogen in a manner similar to that given by Villalobos and 

(5) 

Nuss . A Linde 5A molecular sieve colume was used with argon as the carrier 
gas. 


(4) Experimental Data 


Gas samples were taken from three longitudinal stations along the com- 
bustion chamber. For each of these stations a series of helium injection points 
along the chamber centerline v/ere investigated. 

A typical helium concentration profile is presented in Figure 2. The 
experimental data show considerable scatter, a point which will be discussed later 
in the report. A "best fit" symmetrical curve was drawn by hand through the 

data points to represent the helium concentration distribution. The mean 

. . • — T 

square dxspersxon radxus r of the spread of the helxum tracer gas at a sample 
station was obtci“;;ad from the helium concentration distribution according to 
the relation 


m,, r dr 
He 


( 1 ) 


m„ rdr 
He 


that is , by taking the zero moment and the second moment of the helium con- 
centration distribution at the sampie station. For each of the helium 
concentration profiles, a value of was calculated using Equation 1. The 
values of as a function of the distance between the helium injection point 
and the sample station are plotted in Figure 3 for the three sample stations 
investigated. These results were used later to determine the intensity of 
turbulence and the Lagrangian correlation coefficient. 

The experimental results did not show any significant variation in the 
spread of the helium concentration with the circumferential angle of the 
sample probe. This shows that the assumption of the circumferential homo- 
geneity of turbulence in the combustion chamber is satisfactory. 
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III. ANALYST? OF THL’ EXPERIMENTAL RESULTS 


(1) Steady State Assunption 

As was tTientioned earlier in this report, the data exhibits considerable 
scatter. An accuracy analysis has been perfonned ard it indicated that the 
scatter in the data should be an order of mapnitude less than was observed. 
Also, all of the steps in the sample patherinp and analysis procedure have been 
thorouphly checked. It has been concluded that the scatter in the data re- 
present a true measure of the variation present in the combustion chamber. 

Since the sample time was about one second in duration this would indicate 
some type of very low frequency instabilitv. It might even mean th«t the 
steady state burning pattern in the cham.ber is in some way dependent on the 
ignition history. 

A basic assumption of the analysis herein is that a time average of the 
diffusion of the tracer gas can be obtained by averaging measurements from 
several different runs. This is equivalent to assumir,' steady state turbulent 
flcv as opposed to non-steady flow. 


(2 ) Equations of the Turbulent Diffusion 

In analyzing turbulent diffusion, G. T. Tavlor used the Lar'»'anginn 
approach by considering the path of a marked fluid particle during its motion 
through the flov; field, and develop 3d an equation relating the disolacerent 
of the particle to the turbulent velocity. This formulation was based on the 
assumption that the turbulence was isotropic and homoreneous throughout the 
flow field. 

For the present work, Taylor’s turbulent diffusion theorv was employed 
to analyze the turbulent diffusion in the combustion chamber. However, the 
original assumption was extended to consider that the turbulence field was in- 
homogeneous in the axial direction. Based on thl assumption, the equation 

■1 

for the mean square dispersion radius t in the cylindrical coordinates can 
be written as 


2 

r 


4 


C 


dt' 



R, (T ) 
1-1 


(f) 



-dt” 


( 2 ) 


4 


'■''.'"n t' '’rd twr th<* 'Itrld onrtTrlr and 

X = t’’ - Th'' rapranriar. ronrrlat^'^n rrn^ '■icinn-*- R. (7' *ri t»ic e-’i 

Lj 

is ’ a': 

n /_N - ~ • > v'(t") 

^'i; •• ~~_.zr:- 

(3) 


^ (f) S v’^ ft") 


In order* to solve Equation 2, it is neoess'*r\' that the Larnanrian cor- 
relation coefficient be ropno«;pn'‘ed bv a suitable function. ?’ean”“'OTr.cnts of 
turbulent flews by nany reroa^’chers have shown tha•^ the La^rnanoian cor- 
relation coefficient can for many cituati''ns he rer^osented by an exponential 
function. In particular, Taylor also ured an c>:vcncntial function in bis 
formulation. For the present v.'ork, the Layrargian correlation coo: f icient was 
approximated by the form 


R. (t) = e 


-r.T 


(4) 


where c> is a constant. 

The variables rf and^, v’ * in Equation 2 are functions r " the dh-rersion 
tine. Since the "low iii the corhustion chamber is con'^idered to be* stead*', the 
dispersion tine, t, required for a fluid particle to travel an axial distan<'e 
from the point of release, Z^., to anv point of *ntcrest, Z is related to the 
axial coordinate by the expression. 


t = 



(5) 


where u is the mean axial flov; velocity. Therefore the variables n''*ntioned 
above are also functions of Z. 

The mean axial flow velceitv was not ncasu’^ed in the experiment but vfas 

calculated analytically by the nu'*crical method of the CIC! conpute* pmorran 

(4) 

developed by Rocketdyne , The calculated results of the mean axial flov; 
velocity are presented in Fif'U'''e 4. 


(3) Detorm’nation the Turbulence Fararetors 

Usinm the experimental val’ies for the mean snuare disromnion radii 
(Figure 3), Equation 2 can nov; be solved ^or the root mean scniare (rr?s) 
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600 

(182 



Axiil Distance - in (cm) 





turbulent velocity using a numerical method. The calculation procedur- it* 
given below. 

With the mean square dispersion radius as a known variable and assuming 
a value for a In the exponential expression for the Lagranglan correlation 
coefficient, a first approximation of the rms turbulent velocity as a func- 
tion of the z coordinate was calculated from Equation 2. Then the calcula- 
tion procedure for the equation was reversed. Considering the rms turbulent 
velocity as a known function and using the same value for a as before, tl>e 
equation was solved for T2. The calculated and experimental values of r^ were 
then compared. This calculation was repeated by adjusting the values of the 
rms turbulent velocity until satisfactory agreement between the calculated and 
experimental values of r- was achieved. The above calculations were performed 
for a equal to 100, 4,000, 10,000 and 20,000 sec The comparison of the cal- 
culated and experimental values of r*- Is shovm in rigure 5. 

From the results shown in Figure 5, it can be seen that the curves for 
the three larger values of a investigated, i.e. a equals either 4,000, 10,000 
or 20,000 sec give agreement within the accuracy of the experimental data. 

On the other hand, a rms v'-'oeity curve could not he found which would fit the 
experimental data when a “ 100 sec For comparative purposes herein values 
of a “ 4,000 sec ^ and 10,000 sec ^ will he used. The rms turbulent velocity 
for the combustion chamber using these two values of a are shown in Figure 6. 

For an isotropic turbulence field, the intensity of turbulence T is 
defined as 

T - \ v'" /u (6) 

and the intensities of turbulence for these two values of a are shown on Figure 
7. The results indicate that, for n “ 10,000 sec the intensity of turbulence 
increases from the propellant injector to a maximum value of about 14% at about 
7 inches from the injector, and decreases to about 11% near the nozzle. 

The results of the present work were compared with the results of the 

(2) (3) 

previous investigations using LOX-heptane and GOX-CH^ as propellants. The 
three combustion chambers were identical except for the propellant combination 
and injector used. Figure 8 shows a comparison of the intensity of turbulence 
for these cases. Since the intensity of turbulence is dependent on the value of 
a, any comparison should be made using the same value of a. Since a value of 
a “ 4,000 sec ^ gave satisfactory agreement with experimental data for each of 
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the three configurations this value was ufifd for cooq;>arative purposes in 
Figure 8. As can be seen the intensity of tui oulencc Is greatest for the 
LOX-heptane configuration and is leacr for the GOX-GH2 configuration except 
in one small region near the injector. 

Figure 9 shows the mean square dispersion of a gas particle which is 
released at the injector face and allowed to disperse for the complete length 
of the combustion chamber for each of the three configurations. As would be 
expected the LOX-heptane shows the greatest dispersion while the GOX-GH^ has 
the least dispersion. 
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Fig’jre 3 Comparison of "ntensity of Turbulence 
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CONCLUSIONS 


The diffusion of a tracer gas In a LOX-GTI^ c’oeket «.on^ustlon c'lanber 
has been measured, and the data analyzed using Taylor's turbulent diffusion 
theory to determine the intensity of turbulence and the Lagrangian correla- 
tion coefficient. The results show that the intensity of turbulence 
reaches a maxinium about 7” downstream cf the injector and decreases only 
sligh'^ly to about 122 just upstream of the nozzle. The Lagrangian coirela- 
tion coefficient can be represented by an exponential function of the form 

— OIT *•! 

e w!iere t ■ 10,000 sec 

Among the three rocket combustion chambers investigated in the series 
of experiments at Tulane, the LOX-heptane rocket using a lihe-on-like 
i'<jector produces the highest intensity of turbulence. For the LOX-CH2 and 
GOX-GH2 rockets using coaxia’ injectors, the maxlTnum intensities of tur- 
bulence show only a small difference: 152 for the COX-GU^ rocket and 14% 

for the LOX-GH2 rocket, although the overall turbulence in the COX-GH2 
rocket was considerably less than in the LOX-GH2 rocket. Furthermore, the 
maximum intensity of turbulence oc:curred farther downstream from the injec- 
tor in the LOX-GH2 rocket than that in the GOX-Gll^ rocket. 

Scattering of the helium concentration data persistently occurred 
throughout the experiments using GOX-GII2 and LOX-GH2 propellants, despite 
efforts to improve the accuracy of the data. The scattering is believed to be 
due to an instability of the flow in the combustion chambers. 
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